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Abstract 

Thin liquid films are nanoscopic elements of foams, emulsions and suspensions, and form a paradigm 
for nanochannel transport that eventually test the limits of hydrodynamic descriptions. Here we use 
classical dynamical systems characteristics to study the complex interplay of thermal convection, 
interface and gravitational forces which yields turbulent mixing and transport: Lyapunov exponents 
and entropies. We induce a stable two eddy convection in an extremely thin liquid film by applying a 
temperature gradient. Experimentally, we determine the small-scale dynamics using the motion and 
deformation of spots of equal size/equal color, we dubbed that technique “color imaging velocimetry”. 
The large-scale dynamics is captured by encoding the left/right motion of the liquid directed to the 
left or right of the separatrix between the two rolls. This way, we characterize chaos of course mixing 
in this peculiar fluid geometry of a thin, free-standing liquid film. 
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I. INTRODUCTION 


Chaotic or turbulent mixing is essential for many industrial processes, so a profound un¬ 
derstanding is good for applications. However, despite the basic mechanisms for mixing in 
dynamical systems [8, 19] are well understood, mixing characterization in experiments may be 
hard due to the real-world applications. Here, the restrictions are hnite time, hnite length and 
complex geometries. The system under consideration is an aqueous thin him, driven thermally. 
The him is almost two-dimensional and can show a thick phase (several fim) and a thin one 
(~ 10 —50 nm), both immiscible due to the forces separating them, similar to bubbles (the thin 
phase) in a 3D liquid phase. The thinning speed of the him, i.e. the transition of the whole him 
to the thin phase depends primarily on the mixing of thin and thick phase. Without mixing, 
such a him typically undergoes thinning within several hours, with thermal driving, a how is 
established that mixes the phases and leads to a thinning in a time of the order of seconds. 
Consequently, it is of high interest to understand this mixing and eventually quantify it. In 
the above mentioned article [30] the basic physics are explained. Here, we give a quantitative 
analysis of the mixing properties of the experiment in terms of Lyapunov exponents and en¬ 
tropies. We hrst discuss briehy the properties of thin liquid hlms, followed by an explanation 
of our approach to the characterization of mixing for this highly sophisticated system. 

Thin him dynamics is governed by gravitational, capillary and interfacial forces, where the 
latter are specihed in the disjoining pressure. Combining long- and short-range molecular forces: 
electrostatic, van der Waals (vdW), and steric forces [7, 17], the disjoining pressure depends 
strongly on the distance between the interacting surfaces. Whereas hlms on substrates are 
established in industry and research [21], freestanding thin liquid hlms still provide a challenge 
in experiments and theory alike. Consequently, the study of foam hlms is central to current 
scientihc activities, e.g. [6, 10, 15, 20, 26, 32]. We contribute by presenting a way to quantify 
mixing of vertically oriented, freestanding, thermally forced, noneguilibrium foam hlms. 

As a result of the aforementioned force balance two stable equilibria may occur, depending 
on the chemical composition of bulk solution and chosen surface active agents (surfactants): 
Common Black Films (CBF) with a thickness of more than 10 nm are formed when electrostatic 
interactions balance the dominant van der Waals force, and, of course, gravity and capillarity 
[7, 27]; Newton Black Films (NBF) are stable with a thickness of less than 10nm, due to 
repulsive short range steric forces [12, 13]. In this study we will focus on hlms in their transient 
phase before reaching equilibrium with a typical thickness in the range of 0.1 — 1 {am. The ehect 
of additional forces has been studied by a several authors, mostly for micrometer-thick systems 
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[23, 33], 

Quantifying the degree of complexity of an evolving system is an ubiquitous problem in 
natural science [4], We will focus on two ways of measuring dynamical complexity: the metric 
or Kolmogorov-Sinai entropy kxsi which measures the rate of information production as a fluid 
particle evolves along a pathline, and the Lyapunov exponents (LEs), which give the rates at 
which nearby fluid paths diverge. The principle concept of the KS-entropy is very natural, as 
the information contained in the time evolution is a characteristics of the underlying dynamics, 
cf. Brudno’s theorem [5]. 

From data, one can obtain an estimate by studying the symbolic dynamics acquired by 
assigning different symbols to different “cells” of a hnite partition of the phase space. The 
probability distribution of realized sequences (words) is a signature dynamical evolution. The 
average information gain is obtained by comparing sequences of length m and m+1, in the limit 
of large m: letting the length of the words, m, to inhnity and the partition diameter to zero, 
one obtains the KS-entropy, which is often used as a “measure of complexity” of a system. 

Lyapunov exponents A characterize the exponential divergence of nearby trajectories, typical 
for chaotic systems [14, 18]. One can alternatively state that LEs characterize the sensitivity 
of a system to initial conditions Xj. They are related to the KS-entropy hxs by the Pesin 
formula: 

n 

h = '^\i > hxs ■ (1) 

i=l 

Positive Lyapunov exponents indicate that solutions diverge exponentially on average, negative 
ones indicate convergence. They are computed from time series using embedding techniques 
[14], in our case, we can use the spatial information directly and use instead stretching and 
folding of a fluid area to estimate the local dynamical characteristics directly from the experi¬ 
ment. 


II. THE EXPERIMENT 

In this section, we describe in detail the experiment, the flow structure which is crucial 
for the analysis. Then, we explain the data analysis we applied, hrst for the LEs, then for 
the entropies. In principle, we follow complementary approaches: the LEs are computed from 
microscopic information, whereas the entropies are determined from the large-scale circulation. 
Ideally, both quantities should coincide, to be tested by Pesin’s formula. We will show the 
results below. 
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A. Setup 


The experimental setup consists of a vertical rectangular aluminium frame with rounded 
corners, 45 x 20 mm, enclosed by an atmosphere-preserving cell with a glass window for video 
recording, cf. Fig. 1. Thermal forcing is effected by inserting a cooled copper needle (radius 
1mm) at the him center (T = —169°C), the needle enters the cell through a htting hole. 
Ambient temperature was constant at 20 °C, the corresponding Rayleigh number Ra ~ 10®, 
such that the how is clearly turbulent. Please note that for very thin hlms it is not clear how 
viscosity changes with the him thickness, since surface forces may matter. We are not aware of 
any investigations in that direction, and will not discuss this further. The given number denotes 
an order of magnitude, such that only dramatic changes in viscosity are relevant, e.g. if the him 
has a transition to thin black him and changes structure, too. The temperature across the him 
(z-direction, cf. Fig. 1) is approximately constant, and the Marangoni number Ma ~ 0 . The 
solution from which the liquid him was drawn consists of the surfactant n-dodecyl-/9-maltoside 
{/ 3 —C 12 G 2 ), prepared with hltered deionized water and stabilized with 25 %voi glycerin [24]. 
The thin liquid him is illuminated with a dihuse broad spectrum light source and its rehection 
is captured by a high speed camera. Our vertically oriented foam him is produced initially 
thick (500 - 5000 nm) by pulling it with a glass rod from the reservoir. Quickly, a wedge-like 
prohle develops with BF in a small region, with a sharp horizontal boundary towards the thick 
him below, cf. Fig. 2. 

The interference of incident and rehected light yields a striped pattern, which can be used 
to infer the him thickness. Each color cycle (red —)■ blue) corresponds to multiples n of the 
smallest negative interference condition (2 ■ n -|- l)Xr] = 4hcos0, where the refraction index, 
7], is assumed to be temperature-independent; 0 is the angle of incidence [2]. The velocity is 
measured by color imaging velocimetry (CIV) [29] with u ~ 0.02 ™/s. 

For ideal 2D hows -turbulent or not- surface forces are neglected. This is reasonable for 
thick hlms. In our case, a full description involves all occurring forces [30]. For the Rayleigh 
number we are working with, Ra ~ 10®, other forces can be neglected to hrst order. Surfactants 
stabilize the him, which is conhrmed by our experiment, even if the him is heated (giving rise 
to thermal huctuations) and if the entire him is black, and of the thickness of a few nanometers 
(giving rise to sensitivity to other molecular forces). The turbulence generated by the cold 
needle is observable with the naked eye, a snapshot shown in Fig. 2. 
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Figure 1: Setup of the experiment, (a) schematic sketch of the free-standing thin him, (b) 
light source, (c) camera, (d) cooling rod with liquid nitrogen reservoir. 


B. Flow Description 

This publication presents experimental work and subsequent data analysis, thus we do not 
discuss the equations of motion of the him, but refer to recent work and reviews for the suiting 
equations in our situation [9, 16, 17]. However, we now discuss the observations in order to give 
the reader an impression about the how we consider in terms of dynamical systems properties. 

We use a point thermal force to drive our system, such that we have two rolls, one on each 
side of the cooling rod, cf. Fig. 2. The cold rod drives a stable two-eddy convection with a jet 
of huid in the center as the driving mechanism and separation at the same time. The dehection 
of the jet at the lower frame border is sequential and non-deterministic. For the analysis of 
the large-scale how, the alternation between the left and the right eddy/vortex was visually 
tracked, recorded during a long time interval and converted into a binary time series as detailed 
in section IIC. Further information on this technique is available in our previous works [28-30]. 

In an abstract way, our how can be described by geometry (boundary conditions) and four 
hxed points: the how is bounded by the enclosing rectangle, which in the upper region is 
given by the black him region. Boundary conditions are complicated, in a hrst step they are 
assumed no-slip. Then, we have the centers of the convection rolls, which are elliptic hxed 
points There are two hyperbolic points at the bottom and directly below the lowest part of 
the ice surrounding the cooling needle and a separatrix running top-down (cf. Fig. 3), whose 
connection form a separatrix. This description holds on short time scales, for long times, we 
must consider that the whole system is in a transient state which, however evolves on much 
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Figure 2: Snapshot of a turbulent convecting 
thin him. An image of the non-moving him 
has been subtracted, the frame and frozen 
center region including the cooling tip are not 
shown. On the sides, the typical meandering of 
rising thin/sinking thick him is observed. On 
the top, a layer of black him has formed which 
is transported downwards. Thickness can be 
read from the colors. Left: The convection is 
not yet fully developed, such that “corner” 
vortices are observed; they disappear after 
transients. Right: Turbulent convection has 
fully developed and thin him is transported 
and mixed quickly by the convection. 

larger time scales, such that our considerations are reasonable. A look at corresponding video 
material will clarify this characterization. 

We want to characterize the how, based on measurements. Key characteristics for mixing 
and how in general are stretching or folding of the liquid hlaments. We consider respective 
orbits of huid elements advected by the two convection rolls, which are stably positioned below 
the cooling rod. One key question regards the mixing inside the rolls, the other one mixing 
between the rolls. In contrast to other systems, as e.g. the blinking vortex [19], the vortices 
are maintained and mixing happens across the separatrix due to small diherences in the middle 
downward how. Conceptually, we can use the basic ideas of twist maps, as one common example 
for reduce dynamical systems showing mixing, he. we study the left-right transport, i.e. the 
mixing between the two rolls. 

Below, in Sec.11C, we explain how typical stretching rates are extracted from the experiment, 
as described in more detail in [29]. We recall the procedure briehy here, in order to be complete. 
Then, we step to an analysis in terms of symbolic dynamics in that we analyze the series of 
liquid transported left and right, i.e., we reduce how to the huid exchange between the two rolls. 
This way we have a local measure by the stretching and a more global one by the transport 
across the separatrix. 
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hyperbolic point 


Figure 3: Orbit schematic for the presented experiment (only the right half of the convection 
zone is shown in detail). The light gray outer line represents the shape of the two convection 
zones in the experiment. The gray traces represent orbits moving away or towards the center 
of convection. The elliptic orbits are marked with a light blue circle around the stationary 
elliptic point (blue cross) in the center of the convection zones. Green marks the limit 

between converging and diverging orbits. 


C. Data Analysis 

The captured video data is post-processed to enhance colors and contrast. To analyze the 
behaviour of domains of the same thickness the video is converted into a binary image where 
white corresponds to a single color respectively thickness. As the spectrum repeats continually 
this technique is only valid if the overall thickness deviation is smaller than a full period of 
the smallest wavelength. Subsequently in each frame all clusters of the same thickness are 
numbered and consecutively linked through the following frames (Fig. 4). This enables us to 
track the volume, velocity, deformation rate and angular velocity of the moving fluid. For each 
cluster the velocity is calculated using the shift of its center of mass per frame. The deformation 
rate is calculated by determining the change of the scale of the principal components per frame. 
Similarly the angular velocity is given by the rotation of the principal component. Averaging 
over all frames then delivers the spatial characteristics of the flow held generated by the cooling 
tip. The cluster hnding algorithm is operating with a linear backwards memory of variable 
depth. However for now a tracking of one frame backwards is sufficient to maintain connectivity 
of each cluster through all frames. The memory needs to be limited as merging or dividing 
clusters would be considered connected, thereby distorting the velocity and deformation rate. 

The intermittency of the mixing him is captured by the probability P{Ax, At): we calculate 
the distance Ax for hxed At (on logarithmic scale). The result is plotted in Fig. 5. We observe 
basically ballistic transport for small At, as could be expected because the more chaotic small 
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Figure 4: example of the cluster identification for one frame, including center of mass and 

selected principle components (not to scale). 
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Figure 5: Probability P{Ax, At) for different time scales At(unit: frames, framerate: 100 
frames per second). The transport is ballistic over one decade. Larger times are not shown 

due to the insufficient number of data. 

scales are not resolved, and the large scale transitions are not found by our cluster analysis. 

In the generated velocity field (Fig. 6) the fluid is stretched in regions of high velocity and 
compressed when it enters areas of lower velocity. Due to the shear which is present between 
layers of different velocity folding happens. These two processes are the main aspects of mixing 
in a two dimensional fluid. Diffusion processes can be neglected as the Reynolds number is of 
the order of 10^. 

The overall mixing can be characterized within each vortex as an averaging of the stretching 
rate over the vortex area. The global mixing between the vortices is characterized by the 
probability of a fluid element to cross the separatrix. This procedure can be seen in the view 
of effective diffusion, because the averaged turbulent velocity field is a kind of random walker 
with possibly anomalous and space dependent transition rates, which when averaged yield the 
diffusion coefficient. The transition from one cell to the other is the minimal setup, as discussed 
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Figure 6: averaged velocity field of the forced advection. The cooling rod is positioned at 

approximately prod{x,y) = [2mm, 10mm]. 


in [1], 

We use a naive approach to analyze our film: take a spot of material of size Ax, compute 
its time derivative approximatively from the time evolution as finite difference and use then 
averaging to determine the macroscopic properties. This is in contrast to studies with particles, 
where the velocity of relative dispersion is calculated, and only possible because we have already 
a field at hand. From the clusters of the same thickness, we obtain the eigenvalues e\ 2 along 
the principal axes; they correspond directly to the size of the spot numbered i. Averaging the 
e® over time and filtering yields an estimate for the fields ei^ 2 {x, y)- Now we compare that with 
the definition of the diffusion coefficient: 

D = lim + _ 

At-s>0 At 


In our situation the limit cannot be reached, due to sampling and consequently we should use 
methods like the FSLE [3], which are ongoing work. Here we show the results for the finite-size 
spots with At = As fixed to the minimal sampling time. We make of course an error in mixing 
different spot sizes, which we counteract by choosing homogeneous spots of similar size. 

The space-dependent diffusion is then estimated by 

_ {x{t + Ast) - x{tf) 

^est A 

where averaging over many times (~ 10®) is performed. Of course this is a very rough approach, 
but we will see that in a certain time range we obtain reasonable results, and can compare this 
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with the enhancement of thinning. This is, however not the full story: for mixing it can be 
desirable to have faster than normal diffusion, the characterization of a real process involves 
finite intervals At, Ax, and the degree of anomaly is given by the scaling At ~ Ax^, with a = 2 
for normal diffusion, since we prescribe At and determine Ax accordingly, the full statistical 
characterization is given by the probability function P{Ax,At), cf. Fig. 5. 

So, we characterize mixing within the two rolls by an effective diffusion approach and the 
mixing between the two rolls by statistics of the crossings of the separatrix in the middle of the 
setup. 

The fluid itself offers no contrast to track the motion of hlaments. However, the above de¬ 
scribed reflection imaging transforms the relative thickness difference of convected hlaments into 
a thickness map whose evolution can be followed and analyzed. Compared to other techniques 
like seeding beads or injecting ink, which are unpractical due to the resulting perturbation 
of the thin liquid layer, the thickness tracking is non invasive and precise. A drawback (not 
important here) is that only a relative instead of an absolute thickness map is available and 
hne layering of hlaments is subject to dihusion, thus limiting the resolution. 

In this study we investigate as well on the center stream in between the two stable vortices 
which is dehected either to the right or to the left at the bottom of the frame. At the cooling 
rod huid from both vortices is formed into a center stream. As the convection is bound by the 
frame the center jet is dehected sequentially to the left or to the right at the bottom hyperbolic 
point. This gives a relative measure of the amount of huid mixed m,acroscopically between the 
left and right eddy. The binary series of transports can then be used to estimate the efficiency 
of this mechanism by calculating the entropy and stochastic properties. 


III. METHODS AND RESULTS 

In this section, we present methods and results based on the data obtained as described 
above. A hrst subsection treats Lyapunov exponents, estimated by the stretching and folding 
of small huid elements. In this sense we only hnd hnite-size values, which characterize the 
small-scale properties of the how (and not the inhnitesimal ones, as one would need for the true 
Lyapunov-exponents). Within that subsection we briehy recall and comment on the relation of 
measures from mixing theory and dynamical systems. In the second subsection results on the 
estimation of entropies are shown. We estimate the characteristic entropies by the symbolic 
dynamics approach [1, 14] forming words of a certain length, compute their frequencies as 
estimators of the probabilities, and eventually obtain an estimate for the entropy production 
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with increasing word length. At the end of this section we compare both measures using known 
relations. 


A. Lyapunov Exponents 


In dynamical systems theory, Lyapunov exponents (LE) are key quantities for the charac¬ 
terization of systems, an equivalent is the efficiency in mixing theory for fluids. Essentially, in 
a fluid, the stretch of a fluid hlament is important, analogous to the stretching of a phase space 
element in a dynamical system. In our case, both quantities are identical. Furthermore, the 
mixing efficiency is nothing but a normalized version of the Lyapunov exponent for an ergodic 
system, since time and ensemble average need to coincide for the equivalence to hold. 

In the following we describe how to estimate mixing efficiency and LEs A starting from the 
stretch s of a fluid hlament along its trajectory, where 


s = hm —— 

Axo-5>0 Axq 


( 2 ) 


for the huid velocity measured in the Eulerian (i.e. laboratory) frame. Axq denotes the initial 
separation of two “huid points“ at time t = 0, and Axt the corresponding separation at time t. 
To satisfy the limit towards elements Axq of inhnitesimal length, the cluster data was hltered to 
contain only the smallest clusters available. The local stretching rate ^ (cf. Fig.7) is measured 
in a moving reference system (Lagrangian), therefore the velocity of the separating edges of a 
huid hlament are observed. 

The Lyapunov exponents A describe the typical property of chaotic systems, in particular 
2D systems with hyperbolic points: the trajectories of two nearby points diverge exponentially. 
We can relate that to huid motion, for details cf. e.g. [18, 25]. The exponential divergence for 
inhnitesimally small times is expressed as Axt ~ AxqC^^ . The local Lyapunov exponent is then 
given as the long time average of the logarithmic stretching rate (cf. Fig. 8): 


D In s , 1 

eg = —— —)• A; = hm - 
Dt t 


T 

j Csdf' , 


(3) 


where we recognize the diherence to the stretch (eq.2) in averaging the logarithm or the ratio 
of the lengths directly. For ergodic systems, time and (phase) space average are identical and 
we can compute alternatively 

with the notation (for the space average. We can compute the two Lyapunov exponents 
from the contraction (-) and elongation (-|-) of the tracked huid hlaments and obtain as an 
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Stretching rate <ds/dt> 


Figure 7: Semi-logarithmic plot of the spatially-averaged stretching distribution. Larger and 
smaller unity is marked differently for better visibility. The main difference lies in the 
behaviour for contracting values, which can be zero, i.e. the spot observed vanishes. This is 
reminiscent of our data analysis, explained above and is not expected for a perfect 

measurement. 


estimate on the overall mixing properties of the flow: 

(A_>,, = -1.43 , {A+>,, = 0.92 . (5) 

If the flow is conservative |||A_/A+ = 1|||, here we hnd a contraction considerably smaller than 
—1/0.92 ~ —1.087. It is explained by the fact that we do not only measure the stretch, but in 
addition the general volume loss of a spot of a certain color due to him thinning and diffusion 
of the thickness map, as explained above, cf. Fig.4; this contributes to the negative LE. 

The distribution of the positive and negative LE are quite different. This is explained by 
the dynamics: elongation -positive LE- takes place mainly in the rapid how region between the 
two hxed points (the separatrix). This leads to a very pointed distribution with faster than 
exponential decay, cf. Fig. 8, black marked part. In contrast to this contraction happens over 
the whole region of the him with quite diherent and much slower dynamics. Additionally, the 
contraction rate data has a lower signal to noise ratio compared to the elongation rate, as the 
absolute values are smaller and closer to the resolution limit of the setup. 

Before stepping to the estimation of the entropy we want to describe the connection between 
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Logarithmic stretching rate <£ > 

Figure 8: Distribution of the stretching rate e^. The average over time gives the Lyapunov 
exponents = —1.43 , (A+)^^ = 0.92, indicated by the vertical lines. Clearly, the 

distribution of negative and positive LE ist different. 

LE and mixing efficiency, both quantities are closely related: 


{E.) = ^ 


D : D 


-At' 


( 6 ) 


We see that the difference lies in the nhe normalization. It involves the symmetric part of the 
velocity gradient D = ^ Vn + (Vv)^ (equivalent to the stretching tensor), where D : D needs 
to be constant over the pathlines, for this expression to be valid. It is used as a normalization 
to obtain a stretching efficiency smaller or equal to one. Using this normalization, one can 
demonstrate [19] that the upper bound of the efficiency for two-dimensional flows is Eg,max = 
\/2/2 = 0.707 [19]. Since the stretching tensor is not readily available from the cluster data and 
is only used as a normalization in Eq. 6 it was not used here and we refer to the unnormalized 
Lyapunov exponents. 


B. Entropy 

Similar to the Lyapunov exponent the entropy production of a dynamical system is a measure 
of its chaoticity. Here we compute the Kolmogorov-Sinai entropy h^sfor the center jet deflection 
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and compare it to (A_/+) which we obtain from the stretching rates of the fluid hlaments. 
The topological entropy h is the upper limit for average Lyapunov exponent: h > (A) [25] 
and the Pesin formula equates the topological entropy with the sum of all positive Lyapunov 
exponents: h = A+,j > hxs, and sets the limit for hxs [22]. More precisely the relation 
between the topological entropy h and h^s is given via the Variational Principle: h{T) = 
sup {h{fi)Ks ■ h ^ 'Pt{X)), where /i ranges over all T-invariant Borel probability measures on 
X. Thus h is an upper limit for the KS-entropy. A direct relation between Lyapunov exponents 
and entropy is available via the information dimension 


Di 


h{fi) 




(7) 


of an ergodic invariant measure /i of a smooth invertible map with Lyapunov exponents 
(V)., > 0 > (A_)., [18, 31|. (A Lyapunov exponents calculated from stretching 

(+) and contraction (-) rates. 

Comparing Eq 7 to the Kaplan-Yorke conjecture [11]: Dl = 1 + (A+)^^/| (A_)^J , with 
Lyapunov dimension Dl, suggests that = (A+)^^ in the case of the natural measure of a 
2D smooth invertible map. 

Let us now compute the values for the entropy to the degree our data allow with respect to 
length of the time series and accuracy. First, we need to dehne a kind of “alphabet” to identify 
words. Our procedure follows essentially [1] and for more details we refer to that publication. 
The jet is deflected either into the left or into the right vortex at the bottom of the frame. We 
analyze the deflection pattern to check for deterministic components, which would lower the 
mixing efficiency of this fluid transport. 

The sequence of left (treated as 0) and right (treated as 1) transports is a binary time 
series from which one can associate a word of length n, out of a hnite alphabet: hh^ = 
(S'fc, S'fc+i,..., S'fc+n-i). The block entropies Hn are then calculated from the word probability 
distributions P{W^)-. 

Hn = -J2 7"(W"”) In P{W^) (8) 

where represents the set of all possible words of length n. The entropy per unit time is 
dehned as: 


hn — Hn+1 Hn 
hxs = lim hn 


(9) 


In the presented case the practical limit of hxs is given by the hnite series of events thus 
truncating the number of possible words and the possible entropy gain by increasing the word 
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length, hn can be interpreted as the rate of information production, which for a hnite signal de¬ 
creases when the combined length of all possible words hh" becomes larger than the signal itself: 

> ^s- The mixing efficiency is maximal when the jet deflection is equally distributed for 
left and right. Therefore, we use a set of equal left/right probability as a benchmark for the 
measurement data . As expected, the entropy hn remains constant with increasing word length 
until the hnite size limit of the data set is reached, and longer sets hit this limit at higher word 
lengths, cf. Fig. 9a. 

The entropy of the different distributions below this limit is identical, in particular for 
the data with the same length as the measurement set. Therefore, the deviation in hn of 
the measurement data when compared to the equally dstributed data sets cannot simply be 
omitted as a hnite size ehect. hn of the measurement data is approximately 15-20% below the 
possible maximum given by the random sequences, which indicates a non-uniform process or 
deterministic behaviour, cf. Fig. 9b. 

One might assume that the sequence is tilted towards one side and that this preference 
causes the how pattern to be more predictive. To check for this property we compared the 
data set to skewed random distributions which have a lower entropy production due to the 
higher predictability of the signal. In Fig. 10 random sets with an uneven distribution up to a 
probability of P = 0.9 for one direction of the how are compared. The entropy production of a 
77/23-skewed distribution is comparable to the measurement data. In this case the dehection of 
the center jet would be directed towards one side approximately 8 out of 10 times. However, the 
distribution of left and right transports in the measurement data is even which shows that the 
deterministic components of the signal is not due to a simple asymmetry in the experimental 
setup. 

conditional probabilities Another option to check for deterministic components is to look 
for recurrent sequences in the signal. Therefore, we calculated the conditional probabilities 
P{W^\W'^) with n = m = 1..4 to look for recurring patterns with a memory of up to 4. This 
limit applies due to the hnite data set size, therefore the total word length Iw = n + m ioi 
our measurement is restricted to 8. For each word length the conditional probabilities of all 
possible combinations were calculated giving a matrix of size n x n. Again we compare to a 
random signal of the same size as the measurement data set. For a random data set of inhnite 
length one expects an even distribution of conditional probabilities. However, as we compare 
hnite size data sets, the conditional probabilities become nonuniform at higher word lengths as 
not all combinations are represented equally, cf. Fig. 11b. The hgure displays the deviation 
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Figure 9: Block entropy Hn and entropy hn with increasing word length of the measurement 
set (set length 2 ■ 10^), compared to random distributions of variable set length. 
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Figure 10: Block entropy Hn and entropy hn of the measurement data compared to random 
distributions with varying skewness. Data set length is constant. 


from a perfect random distribution, which would show as a 50% gray tile. 

The conditional probabilities of the data set deviate substantially from the even distribution 
the random data provides, cf. Fig. 11a. The jet is more likely to alternate between left and 
right transports which is evident by the higher conditional probabilities of alternating binary 
combinations, ie. P(fF”'|1010), P(fF”|0101). Some combinations with two or more consecutive 
transfers in the same direction show a higher than average probability, but these are always 
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combined with alternating sequences. Uniform combinations, ie. P(1U'^|0000), P(1U'^|1111), 
are extremely rare events. 

The diminished entropy production with increasing word length is caused a preference for 
alternating transfers which are evidenced by nonuniform conditional probabilities. The signal is 
more predictable and most common combinations are already covered with lower length words 
thus the additional words carry less new information. 



(a) Measurement data (b) Random distribution 


Figure 11: Conditional probabilities for measurement data (a) compared to a random 
distribution (b). Each tile represents the conditional probability P(1U™'|1U”), where lU™' is 
plotted on the x-axis and IF” is plotted on the y-axis. Word length = 1..4. The grey scale 
for each plot is centered at the expected probability, i.e. for Iw = 2: P(1F™|1F”) = 0.25 for all 
possible combinations. From this value the scale covers 50% deviation, i.e. Iw = 2: black: 
P{W^\W^) = 0.125, white: P{W^\W^) = 0.375. 


Eventually, we compare the results from entropy and LE estimation. For the LE we obtained 
A+ = 0.92, A_ = —1.43. For the Shannon entropy we can read off the entropy gain from Fig. 10 
(red line and markers) as Hks = 0.82. Using the Kaplan-Yorke conjecture we should hnd 
a coincidence h ~ (A+) with a Kaplan-Yorke dimension of 1.92. Given the error sources in 
the LE estimation and the relatively short time series, we consider this result as a very good 
coincidence. 
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IV. DISCUSSION AND SUMMARY 


We presented a relatively simple experiment exhibiting complex dynamics, where turbulent 
mixing reaches out for the nanoscale, at least in one dimension. This renders the flow two- 
dimensional, however with additional forces acting between the surfaces: disjoining pressure 
and capillary pressure. 

With respect to turbulence, we observe that the relatively low-i?a convection generates 
a weakly turbulent flow with two prominent rolls, the turbulence mainly can be observed 
inside the rolls on shorter time scales, whereas transport between the rolls is on a slow scale 
and shows signs of chaos. To characterize the flow, we take advantage of our measurement 
capacities, and the color imaging velocimetry (CIV) technique. In addition to the flow held, 
concluded from Lagrangian trajectories, we obtain local deformation rates which are used to 
estimate the Lyapunov exponents. With respect to mixing we clearly observe the typical chaotic 
hlamentation. As mentioned in the context of weak turbulence, global mixing on the large scale 
happens between the left and right convection roll. 

We focused on two different methods to quantify the dynamics of the how: Lyapunov expo¬ 
nents and entropies. Whereas the LE are calculated using small-scale dynamics using CIV, i.e. 
tracking of small huid spots, entropies have been calculated using the left-right transport across 
the separatrix. Qualitatively, these two observations are related by the two “spots” generating 
the dynamics: the hyperbolic hxed points. They determine the microscopic properties (LE) 
and as well the transport across the separatrix. Thus, based on dynamical systems principles 
one expects that both quantities coincide. Given the two diherent approaches and the very 
diherent spatial scales, this coincidence is a formidable cconhrmation of theory. 

The presented data is gathered from 8 individual runs of the experiment. Entropy production 
and conditional probabilities were calculated for the data set and compared to truly random 
data sets of varying size and skewness. The conditional probabilities show that an alternating 
pattern of left and right transports is preferred which increases the predictability and lowers 
the entropy production of the measurement data with respect to equally distributed data. The 
number of left and right transports in the data is even with a deviation of 2% from the exact 
left/right randomly equal distribution, which can be accredited to inaccuracy of experimentally 
obtained data (rod not positioned at the exact center) and limited data set length. Although 
the entropy production is comparable to skewed random distribution with approximately 
70/30 shift, these are two unlinked deterministic mechanisms which lead to a decrease of 

Due to the shortness of the time series, an entropic analysis has errors beond a certain 
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length of the words formed. In order to have a quantitative measure for this length we compare 
the measurement with numerically determined left/right random sequences without memory. 
As displayed by the random data sets of varying size a larger data set would yield a higher 
accuracy in the entropy production rate and conditional probabilitiyes, but the fundamental 
scaling behaviour remains the same. Thus the statistical evaluation of the measurement data 
is valid. 

Using such kind of numerical validation procedure, we hnd clearly that the asymmetry 
found should be attributed to higher-order memory in the data. This is quite plausible, since 
we are investigating a fluid system, where memory is “built in” for advected particles and so 
for advected fluid particles, too. The calculated KS-entropy serves as an upper limit to the 
Lyapunov exponents, which were calculated from the stretching rates. Experimentally, we used 
color-tracking which focuses on small scales. This area tracking algorithm provides detailed 
information about the flow held and is adaptable to work on any data of deformable clusters 
with high enough contrast with respect to their surrounding. 

One key question regards the consistence of LE and entropy estimate. We found an astonish¬ 
ingly close result for entropy and LE using Kaplan and Yorkes conjecture. Both agreed within 
an error of 10%. The values of 0.92 and -1.43 indicate a relatively moderate large-scale mixing 
which is due to the fact that the setup involves only two main how regions among which one 
has to achieve mixing. 

However, in practical terms, such nanohuidic devices can be well used as free-standing mixers 
with very hexible surface properties. The mixing times are then quite fast, since the device 
dimensions can be miniaturized further. The true advantage lies, in our opinion, in the very 
good ratio of surface and volume of the him: such a him has an enormous surface, in our setup 
the ratio was 1:100000, which is probably only achievable with thin hlms. The free-standing 
property, hnally, allows to avoid problems with solid-liquid interface forces, rather one can focus 
on a chemically favorite design of the equally treated surfaces. 

Generally, we think that with our relatively simple experiment we can touch deep physical 
questions on the microscopic nature of thin him hows and treat at the same time very practice- 
oriented topics. 
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